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Using effective meson Lagrangians we study the interaction of rho mesons in a hot baryon- 
free system. Various mesonic resonances in direct s-channel reactions are investigated employing 
standard self-energy techniques, including new reactions that have up to now not been considered in 
a self-consistent approach at finite temperature. The importance of subthreshold resonances, which 
are readily accounted for through off-shell effects within our framework, is emphasized. Special 
care is taken in reproducing radiative decay widths, as they provide valuable constraints on the 
evaluation of dilepton spectra. In particular, we compare our results for dilepton production rates 
... to earlier calculations based on an incoherent summation of individual processes. 
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D-l' I. INTRODUCTION 
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\Q ' Dilepton measurements in heavy-ion collisions at intermediate Q and high |g,|| bombarding energies have revealed 

a strong excess of pairs as compared to proton-induced reactions. It has become clear that the collection of hadronic 
sources that can successfully account for the measured spectrum in proton-induced reactions fails to reproduce the 
^^ measured yield in nucleus-nucleus collisions. Even after the inclusion of (free) 7r"'"7r~ annihilation during the interaction 

phase of the hadronic fireball the enhancement, especially for invariant dilepton masses Mu below the p-mass, remained 
unexplained B . Thus many theoretical efforts have concentrated on the role of medium effects in pion and rho meson 
/^ propagation. It seems fair to say that there currently exists two main schools of thought that have staked a claim 

fS| onto the theoretical interpretation of the CERN low-mass dilepton measurements. The first one assigns the rho meson 

^^ mass as an order parameter of the chiral symmetry restoration. The CERN data, in this approach, then signals a 
0^ dropping of the rho meson mass mn. A second interpretation relies on the fact that in a strongly interacting medium, 
^^ , the rho meson will have its width greatly increased due to modifications of its pion cloud [p-13| as well as to the 
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f~| ■ direct coupling to baryonic resonances |14|-[l7||. Whether these two different scenarios can be reconciled, or whether 

OhI one of them can be eventually ruled out, is the subject of much current research and debates. In this article we follow 

Jv • the approach that is germane to the second of those theoretical avenues. We will restrict ourselves to a heat bath of 

fli ] pions, kaons and rho mesons, characterized by a finite temperature T. Using the many-body formalism of Ref. p5[ , 

r^ , we will separately investigate each meson channel and calibrate its individual strength through empirical information 

L| ' on both hadronic and electromagnetic branching ratios, thereby introducing some new channels that have not been 

. ^H considered before in a framework like the one at hand. 

^^ The baryon-rich nuclear medium has been shown to strongly affect the vector meson properties, see, e.g., 
Refs. |l5| , [l4| , [12 171^. Then, a self-consistent finite-temperature assessment of the meson contributions is neces- 



sary to quantify this assertion. This is particularly true for CERN-SpS energies: in spite of the fact that the phase 
space is dominated by mesons (with a final pion-to-nucleon ratio of about 5:1), the explanation of the dilepton data 
in both the dropping mass and the in- medium broadening scenario crucially depends on baryonic contributions. At 
RHIC energies, however, it is hard to imagine that baryons could play a major role; since dilepton measurements in 
the PHENIX experiment will be able to address the invariant mass region between u;(782) and (/)(1020) with excellent 
resolution, in-medium properties of the rho in mesonic matter should be most relevant there. Various analyses in this 
direction have already been performed po|- p4| , [15|j2^ ] . The moderate collisional broadening for on-shell rho-mesons 
found in Ref. | 
results were o 



was mainly attributed to resonant scattering via intermediate ai(260) and ^^"1(1270) states. Similar 
Dtained in Ref. |1^ using finite temperature self-energies which include off-shell effects. This is to be 
expected, since the ai(1260) and Ki{121Q) are located above the free pir and pK thresholds, respectively, which 
makes the resonant contribution dominant. In order to properly address the low mass region, additional mesons have 
to be considered. Those larger meson ensembles have in fact been included in calculations where individual rates 



are summed over 21 1. There, the contribution of the uj, for example, appears as a radiative decay channel. In the 
language used in this work, the lj is a sub-threshold state, i.e., in p-K — > u;(782), where, given a typical thermal pion 
energy of 300-400 MeV, the relevant p-mass to form the to would be substantially off-shell, M ~ 400 MeV. 

The objective of this article is twofold: on the one hand, we would like to examine if a consistent off-shell treatment 
of the most important mesonic resonances has a more severe impact on the dilepton production rates than what has 



been estimated so far within different frameworks. At the same time, we can investigate possible interference (or 
collectivity) effects in the coherent summation of the various self-energy contributions. Since many-body calculations 
of this nature can be rather convoluted, we will exhibit an explicit channel decomposition of the considered resonances 
R = U!,hi,ai, Ki,tt' and /i. An upper mass limit of 1.3 GcV enables us to adequately address the phenomenology 
that concerns us in this work. As the simplest, but by no means negligible, finite temperature effect in the pion 
cloud of the rho we will furthermore include the Bose-Einstein enhancement in the p —^ tttt decay width, which is 



not new [E6 20 1. Our paper is organized as follows: the next section introduces the hadronic Lagrangians used in 



this work, followed by the determination of the parameters in our model in sect. III. We then evaluate the resulting 
in- medium rho properties (sect. IV) and their impact on dilepton production rates (sect. M). Sect. VI contains a 
summary and concluding remarks. 

II. INTERACTION LAGRANGIANS 

Our starting point is the model for the p-meson in free space employed previously in Refs. ||l^,|l^,|l^ . Based on the 
standard pmr interaction vertex (isospin structure suppressed) , 

Cpn-^ = gpT,-,, n p^npp , (1) 

(p'^: pion momentum) the bare p-meson of mass mPp is renormalized through the two-pion loop including a once 
subtracted dispersion relation, giving rise to the vacuum self-energy 

t%AM) = / ^ vp^pf G%iM,p) , (2) 

with the vacuum two-pion propagator 



Gl^iM,p) - —^ ; Lo^ip) = ^ml+p^ (3) 

LU^[p) M^ - [2u;^[p)y + IT] 



and vertex functions 



Vp7T-^{p) = v/o 9p7T-R 2p Fp^^{p) (4) 

involving a hadronic (dipole) form factor Fpj,-^ |]l3| (cf. also Eq. (n2|) below). Resumming the two-pion loops in a 
Dyson equation gives the free p propagator 

D°iM) = [M' - (mir - ^%AM)]-' , (5) 

which agrees well with the measured p-wave tttt phase shifts and the pion electromagnetic form factor obtained within 
the vector dominance model (VDM). 

To calculate medium corrections to the p self-energy in a hot meson gas, we will assume that the interactions are 
dominated by s-channel resonance formation. At moderate temperatures relevant for the hadronic gas phase, the 
light pseudoscalar Goldstone bosons P = tt,K are the most abundant species. We can group the various resonances 
in pP collisions in two major categories, namely vector mesons V and axial-vector mesons A. For the latter, a simple 
interaction Lagrangian, compatible with chiral symmetry and electromagnetic current conservation, is given by 

CpPA - GpPA Ap, (g'"- q^p^ - qf^p") p, P , (6) 

although other choices are possible [pSj . pP scattering via intermediate vector mesons V is determined by Wess- 
Zumino anomaly terms, which are of unnatural parity and involve the four dimensional antisymmetric Levi-Civita 
tensor e^'i^'^T. 

Cppv = Gppv e^.ar k^ V q^^ P . (7) 

In both Lagrangians (|g) and (0), p'', qf^ and fc^ denote the four- momenta of the pseudoscalar, rho- and (axial-) vector- 
mesons, respectively. As a third possibility pP scattering can proceed via a pseudoscalar resonance; here we restrict 
ourselves to the process pn -^ 7r'(1300), which can be described by 



l^pPP' - Gppp, P' {k-qVy.-V(lK) P"" P ■ (8) 

In addition to that, we will need a pV A interaction vertex, which is also related to anomaly terms |M. We choose 
the following form, 

LpVA = GpVA tpuar P^V^ P"" fca^" - \ {kpA^f , (9) 

which again satisfies the appropriate conservation laws, p""" = q" p"" — q"" p" is the usual field strength tensor. We 
have explicitly written here the kinetic energy term of the axial vector field where the constant A represents a gauge 
freedom connected with the axial- vector field 1291 . In what follows we take A = 1 ; this choice will be further motivated 



in sect. IV 



An important hint on the importance of resonances in dilepton production is provided by their radiative decay 
width, which constitutes the M^ = (photon) limit of the timelike dilepton regime. In Table | we have collected 
mesonic resonances which are accessible via p-induced excitations and exhibit substantial decay rates into final states 
involving either photons or rho- mesons (or both). We expect these to be the relevant contributions to dilepton 
production. In particular, note that the /i(1285) has a large radiative decay width of 1.65 MeV for /i -^ p7, which 
led us to the consideration of the ppj\ interaction, Eq. (pi) . 

As a first step one has to assign realistic coupling constants in the effective Lagrangians. They are adjusted to the 
experimental branching ratios of the resonances into, if available, both ph and 7ft. (h — tt, K^ p). For reliable estimates 
of the pP decay widths it is important to include the finite width of the p (in particular for sub-threshold states like 
the cj(782) or /i(1285)). This is accomplished by folding the expression for the width at given p-mass M with the 
spectral function A'l{M) — — 2ImD°(Af). For the axial-vector meson resonances in pP scattering the vertex of Eq. 
leads to 

GIpa if (2/p + 1) f MdM 

^■^-"^('^ - ^^ (2/^ + l)(2J^ + l) J -^M^nic^ 

X [i(s - M^ - mlf + M^LOp{q,^f] FpPAiqcm? , (10) 

and from Eq. (M) one obtains for vector resonances 



M" 



,._Glpy IF{2Ip + l) f MdM , 

2m„ 

with qcm being the three- momentum of the decay products in the resonance rest frame, u>p{qcm)'^ = "^-p + Qcm- ^^ 
is an isospin factor, M™°^ = ^/s — mp, and Fppp {R = A, V) are hadronic form factors that reflect the finite size of 
the fields that appear in the effective vertices. We take them to be of dipole form, 

/ 2A2 -I-toI \^ 

Fppp{q.„.) = — , , ^ ^ , (12) 

Y2A^p + [Ujp{qcm) + i^p{qcm)\ J 

normalized to 1 at the resonance mass rap. With the Lagrangian for the pseudoscalar resonance P' — 7r'(1300), 
Eq. (H), one arrives at 

For the pVA vertex, which in our case corresponds to the /i(1285) -^ pp decay, the spectral functions of both outgoing 
p mesons have to be integrated over. One has 

_ Glpf^ nIF{2Ip + l) f A'hdAh ,0.,. . f M,dM, , 

^^^-'"'(^)-^^(2/,,+l)(2J,,+l) J ^^^p(^^i) J -^MM2) 

X 2sgL-F^PP/i(9cm)' , (14) 



where the factor n = 1/2 ensures the proper symmetrization of the two p mesons in the final state. The integration 
Umits are Mf"^ = Vs - 2m^, as before, and M^""" ^ ^ - Mi. 

For each vertex one is left with two unknowns: the coupling constant G and the form factor cutoff Ap, which should 
lie in some reasonable region for the hadronic scales at hand, typically Ap < 1 — 2 GeV. The hadronic decay widths are 
in fact not very sensitive to Ap, since the dominant contributions in the Af- integrals are centered around the p-peak 
in j4"(A/) (or towards the maximal M if mn < nip), where the three-momenta are rather small. More stringent 
constraints on Ap are imposed by the radiative decay widths, since the massless photon can carry away the maximal 
three-momentum. Invoking the phenomenologically well-established (especially for purely mesonic processes) vector 
dominance model (VDM), the photon decay widths follow from the hadronic couplings to vector mesons by simply 
(i) taking the M^ -^ limit, i.e., substituting A^{M) = 2ttS{M'^) for real photons, (ii) supplying the VDM coupling 
constant (e/g)^ ~ 0.052^, and, (iii) omitting the (2Jp -I- 1) isospin degeneracy factor for the final state. This yields for 
both axial-/vector resonances {R — A, V) 

^^^^^ ^ -^ [l) {2In + l){2Jn + l) ^"'"^ ^''^«(«^")' ' ^^^^ 

whereas for the /i — > "fp decay one still has to integrate over one p mass distribution, 

^'^-^"-^^ UJ (2//,+l)(2J/,+l) J ^^Ap{M,)2sq,^FppfMc^) . (16) 

2m„ 

The additional factor of 2 accounts for the two possibilities of attaching the photon to either of the outgoing p's. 
Within the simple version of VDM employed here, the radiative decay of the 7r'(1300) vanishes. 

III. DETERMINATION OF FREE PARAMETERS 

Let us now discuss the individual resonances, based on the interaction vertices formulated above, in more detail. 

The p7rai(1260) coupling constant has been estimated in our framework in Ref. ||l^ as Gprrai = 13.20 GeV~^ 
taking ||2|] Ta^^^j^p = 400 MeV, rua^ = 1230 MeV and assuming a cutoff Ap^rai — 2 GeV. With these parameters 
the radiative decay width of the ai turns out to be 1.23 MeV, somewhat larger than the only available experimental 
information quoting a value of r^j^-^^^ = (0.64 ±0.24) MeV. This can be substantially improved upon with a reduced 
cutoff Ap^a^ = 1 GeV, requiring a slightly higher couphng constant GpTrai — 13.27 GeV~^ to fit the hadronic decay 
width, but resulting in a photon decay width of Fq^^^^ — 0.66 MeV, now in good agreement with the experimental 
value. 

The i4ri(1270), which was also included in Ref. ||l^, is the appropriate resonance in pK scattering. Since there is 
no radiative decay known, we assume (in analogy to the oi cutoff) ApKKi = 1 GeV. Taking Tk^^Kp = 60 MeV, 
■niKi = 1270 MeV yields GpKKi = 9.42 GeV~^, resulting in Tki^Kj — 0.32 MeV, which is qualitatively in hue with 
the naive expectation that in comparison to the ai, the heavier strange quark essentially acts as a spectator in the 
photon decay. 

The /ii(1170) is the isospin-0 pendant to the ai(1260). No quantitative empirical information on its decay properties 
is available. We therefore make the plausible assumption that the major part of its width originates from the only 
observed pn decay channel, i.e., Thi-tTrp — 300 MeV. Using again ApT^hi = 1 GeV gives Gp,rhi = 11.37 GeV~-^ and 
^hi^TT'y — 0.6 MeV, which seems not unreasonable. 

The (x'(782) meson differs from the previously discussed resonances in that it lies significantly below the free irp 
threshold. Hadronic models including the effective four-meson w37r-vertex usually attribute substantial parts of the 
hadronic decay width of about 7.5 MeV to Trp states [|o|. On the other hand, the radiative decay w — > 7r7 should, 
within VDM, entirely proceed through np states. Enforcing the experimentally rather precisely known value of 
^Lu^TT-y = 0.72 MeV, and using again Ap^rt^ = 1 GeV, yields the coupling constant G^^t^ = 25.8 GcV^^, entailing 
F^^^p = 3.5 MeV, which is somewhat on the low side of the typical values |^,^. One could accommodate 
larger values by choosing softer form factors, but more detailed information should be inferred from the dalitz decay 
spectrum cu — > iv^ p,^ p~ , see, e.g., Ref. |p3]. Let us point out here that the Dalitz decay '^u^-nOe+e- is about one order 
of magnitude larger than the direct dilepton decay width F^^g+g- |29|] . 

The /i(1285) is similar to the u;(782) in the sense that it is also a sub-threshold state (about 250 MeV below twice 
the p mass), accompanied by a rather small total decay width of '^ 25 MeV; this supports a tempting interpretation 
as a state with a large pp component, the decay into it being suppressed by the lack of phase space. We will have 
more to say on this later. For the time being, if we attribute the entire decay width into 4 pions of about 7.5 MeV 



to the pp channel, the radiative decay is overestimated by at least a factor of 2. Smaller form factor cutoffs are not 
really efficient in suppressing the radiative decay width, since the latter still involves an integral over one p spectral 
function where dominant contributions arise from masses M2 — rup, where the decay momentum is rather small. A 
reasonable compromise for our purposes appears to fix the radiative decay width approximately at its experimental 
value of T f^_,pa^ ~ 1.65 MeV, which, choosing a form factor cutoff h^ppf^ — 0.8 GeV, results in a hadronic decay 
width Tf^^pp = 3 MeV. 

The total width of the 7r(1300) is not very well known, quoted as 200-600 MeV by the PDG H]. Lacking more 
precise information, we attribute 300 MeV to the Trp channel, which is one of the two observed decay modes. The 
specific form of our interaction Lagrangian, Eq. (0), together with the VDM assumption, does not allow any radiative 
decay. The latter has not been observed so far. 

The coupling constants and cutoff parameters as well as the resulting branching ratios are summarized in Table O. 

IV. IN-MEDIUM p PROPAGATOR 

Having fixed the parameters of the interaction vertices, we are now set to calculate the in- medium p self-energy, Ep, 
and corresponding dilepton production rates. The self-energy is related in a standard way to the forward two-body 
p scattering amplitude off the surrounding thermal mesons. Within the imaginary time (Matsubara) formalism one 
obtains: 

s;a'?o,9;r) = y'^^^[/''K(p))-r''K(p) + qo)] A^^Tfe?) , (17) 

where the isospin averaged p scattering amplitude Mph is integrated over the thermal Bose distribution f'''{ujh{p)) — 
[exp(w;i(p))/T — 1]^^ of the corresponding hadron species h with ujh{p) = \/to^ + jP. The invariant amplitudes are 
evaluated in terms of the s-channel resonance contributions from the previous section. Collisions of the p with the 
pseudoscalar mesons P — 7t,K lead to 



MpPAip, q) = IF G^pA FpPAiqcrafie^ p-q-q'^ p-e) DA,nx{s) [e*^ p-q-q^ p-e*) . 

Then, one obtains for the axial- vector resonances A = hi,ai, Ki : 

KpAiP^I^ = IF GIpa FpPA{q„nf Da{s) «r(P:9) 

v7{P^ 9) = -g'^'^iP ■ qf + q^"^^^^ + [q^p" + qy^P ' 9(1 - -) +P>''<z'(l - -) 

s s s 



(18) 



(19) 



for the vector resonance V = uj: 

Mppvip,q) = IF Glpy Fppviqcm)^^""^'' e^ ^o 1p) DvM^) (e^^^" el k^ qs) , 
M^;^y{p,q) = IF Glpy Fppviqcmf Dv{s) wr(p,9) , 

«r(p, q) = '9'"'{{p ■ qf - p\^) - q'^q'p' + iq'p" + q^^p ■ q - p^'q^ , (20) 

and for the pseudoscalar resonance 7r'(1300): 

Mppp^{p,q) = IFGlpp,FpPp,{qc„if{e -pk-q-e-kp-q) Dp^{s) [e* ■ p k ■ q ~ e* ■ k p ■ q), 

M';pp,{p.q) - IF Glpp, Fppp,{q,^f Dp,{s) v'^^p,q) , 

v^p'^ip, q) = q'^q'^ip ■ qf - [q^p" + 9 V) q^ P ■ q + P^P^'q' ■ (21) 

The amplitude tensors Af ^"^ are obtained by removing the p meson polarization vectors e^, e* from the invariant matrix 
elements and contracting the remaining indices. The intermediate (axial-) vector propagators at four-momentum 
k = {p-\- q) have been taken as [R — A, V) 

n /,N^ {-9K\ + k^kx/s) 

= {-g.x + k^kx/s) Dr{s) . (22) 



Using the same form of propagator for the /i-meson in the pp scattering process raises a problem: the corresponding 
(Born-) amphtude for pp —^ fi —^ pp vanishes identically, which is, in fact, in line with analyses based on chiral 
Lagrangians [ pip^ . However, the latter then suggests, within VDM, the absence of the direct /i -^ p°7 decay, which 
empirically is quite large. Since here we are interested in phcnomenological estimates for dilepton production, we 
decided to circumvent the vanishing ppfi coupling by making use of the gauge freedom for massive vector particles 
provided by the additional — ■|(fc^/f )-term (Stiickelberg term [g8[) in the interaction Lagrangian, Eq. (g). With A = 1 
the /i propagator takes the form 

DJ^is) ^ -g'^^Df.is) , (23) 

which has been used in the actual calculations. We note that when further applying the naive VDM for the two-photon 
decay of the /i, we would obtain a nonzero branching ratio, violating Yang's theorem |3J]. This is, of course, an 
artefact that might be related to our gauge choice (in fact, Brihaye et al. ||3^ pointed out that VDM breaks down at 
the two-photon level). On the other hand, we have found that the /i contributes negligibly to the quantities being 
calculated in this paper, with the parameters delineated as described above. Therefore, the role of the /i, which at 
first seemed promising because of its large radiative decay width, turned out to be unimportant from a pragmatic 
point of view, so that we do not attempt further improvements of its interaction vertex. 

In the medium, the specification of a thermal rest frame breaks Lorentz invariance. As a consequence, the in- 
medium p self-energy tensor is characterized by two independent scalar functions, each depending separately on 
energy and three- momentum (in the vacuum one scalar function depending on invariant mass only is sufficient). This 
is conveniently described in terms of longitudinal and transverse modes of the p propagator |20| : 

^^''('^°' '^ = M2-(mO)2-I]i(go,g) + M^ - (mO)^ -j:T{q,,^ + (mO)2Af2 ^^4) 



with the standard projection operators 



pfuj _ H H _ uu _ pfii' 



, // = or i^ = 

3 



^^'~ <5'^-2^','/i,^G{l,2,3} ^^^^ 



(the spacelike components of p and v are denoted by i and j, respectively). The longitudinal and transverse self- 
energies are defined by the corresponding decomposition of the polarization tensor: 

Y.^;{q^, q) = E^(go, q) P't + Sj(<7o, ?) ^r ; (26) 

they are calculated from Eq. (^7|) as 

^^/^(go, g) = GIpj, if J ^^^^tf^f^p^ iri^piP)) f'^'i^pip) + qo)] F.pniqarrf 

X DR{s)v'i^-^ip,q) (27) 

for R ~ A,V, P' with x — cos 9, — l[p,q). The projected vertex functions are given by 

v\{p,q) = (^l)mi' '"A^P^'i) 

2 2 2 

q t 1 2\2 I ,' \2 2/1 q \ -2 2 2/1 9 \ 

= -T-\s -q - m^) + ^p{p) q 1 )~p q X 1 

4s s s 



VAiP,q) = 2^Pt)^u v'X''ip,q) 



l[hs-q^-ml)'-fq^il-x')il-^)] 
II s 



= q'f{l~x') 



VviP,q) = -^^Pt),,,. v'^''{p,q) 
1 4 



- qis ~q^~ mlf - 2mW - 9^(1 - x^)] (28) 



and similar expressions for the pseudoscalar resonance P'. For the /i(1285) an additional integration over the mass 
distribution of the p from the heat bath has to be performed: 



^L,T 



where M2 = Pq— j? ■ 



'-^A,M J Hg iripo) - /-(po + .0)] F,,,, i,.^r 



(29) 
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FIG. 1. The real and imaginary parts of the polarization-averaged p self-energy (lower and upper panel, respectively). The 
different channels are labelled explicitly and explained in the text. Note that the tttt channel is absent for the sake of clarity. 



In Fig. |l| the real and imaginary parts of the individual spin-averaged self-energy contributions, 



(30) 



(h = TT, K, p) are shown at fixed three-momentum modulus \q\ — 0.3 GeV in the lower and upper panel, respectively. 
Around and above the free mass nip, the strongest absorption is caused by ai(I260) resonance formation, which is 
about as large as the sum of all other channels, shared to roughly equal amounts between iCi(1270), /ii(1170) and 
7r'(1300). The i^i(1270) curve acquires its maximum at lower M than the pion-resonances due to the higher thermal 
energies of the kaons (including their rest mass). In the low-mass region M < 0.6 GeV, the dominant contribution 
is due to the w-meson, which, however, barely leaves any trace in the resonance region. It is also seen that the 
effect of the /i is very small. In the real part of the total self-energy we observe appreciable cancellations, until 



eventually all contributions turn repulsive (the latter feature would of course be modified when accounting for further 
higher resonances). Such cancellations are typical for many-body type calculations as performed here. They are the 
reason that one usually encounters only moderate modifications of the in-medium pole mass. On the other hand, the 
imaginary parts of Sp strictly add up, generating significant broadening. We refrain here from plotting the two-pion 
loop contribution to the p self-energy. It will start exceeding the in-medium corrections past M ~ 0.450 GeV (the 
free p width at M — nip amounts to about 150 MeV). 

As discussed earlier, the existence of a preferred thermal reference frame will break Lorentz invariance. An advantage 
of a theoretical approach like the one at hand is that the transverse and longitudinal parts of the rho self-energy can be 
separately resolved. This is shown on Fig. 0. Even though at the present time one does not have a practical observable 
that is convincingly sensitive to the polarization, one should keep this difference in mind for future applications [gO| . 
From the figure, we see that at finite three-momentum the polarizations differ the most at low invariant masses and 
become undistinguishable at high masses. 
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FIG. 2. The longitudinal and transverse polarization contributions to the imaginary part of the self-energy as arising from 
the sum of meson resonance contributions. 



In addition to the direct p — h interactions we account for the simplest medium effect in the pion cloud of the p by 
including the Bose enhancement in the two-pion bubble. In the Matsubara formalism this amounts two replacing the 
two-pion propagator, Eq. (||), by ||2^ 

1 [l + 2r{uj^{p)-T)] 



G^^{M,p-T) = 



(31) 



oj^{p) NP - {2Lo^{p)f + 11^ 

At a temperature of T = 150 MeV, this generates an additional broadening of the p self-energy starting from the 
two-pion threshold reaching an appreciable maximum of ~ 20 MeV at M ~ 0.6 GeV (on the scale of Fig. 1, upper 
panel) and gradually decreasing beyond. 

To end this section, we plot in Fig. g| the full spin-averaged imaginary part of the p propagator (spectral function). 



Imi?p(M, q\T)^- [lmi?^(Af, g; T) + 2Imi^J(Af, q- T)] 



(32) 



in a thermal meson gas of temperatures T = 120, 150 and 180 MeV as appropriate for the hadronic phase in 
ultrarelativistic heavy-ion collisions. More explicitly, one has 



lmD];-^{M,q-T) 



ImJ:^-^ {M,q;T) 



|M2-(mO)2-I]^'^(M,g;T)P 



(33) 



with the longitudinal and transverse self-energy parts 



'^piTTT + 2_^ '^pa 



^p-rnr + / ^ ^pa i 



(34) 



where the summation is over the mesonic excitation channels a—TTUj, nhi, ttci, tttt', KKi, KKi, pfi, as discussed, 
and 'SpTjj^ now contains the Bose-Einstein factors through Eq. (pl|). We find that the thermal p spectral function 
undergoes a broadening (defined as the full width at half maximum) of about 80 MeV at T = 150 MeV (with little 
three-momentum dependence, see also, e.g., Ref. fl^l), which almost doubles to ^ 155 MeV at T = 180 MeV. Those 
values are a factor of 2 larger than the collisional broadening found in Ref. 122] based on on-shell scattering amplitudes. 
In Ref. [|l8| the p meson self-energy has also been evaluated for on-shell p mesons using the Tph-Qh approximation {i.e., 
the self-energy being proportional to the p-h scattering amplitude and the matter particle density qu). For a pion gas 
of density n^ ~ 1.5 fm^'^ a broadening of 400 MeV has been quoted, which, when rescaling to a density of 0.12 fm~'^ 
(corresponding to thermal equilibrium at T — 150 MeV) gives ^ 30 MeV, again about a factor 2 smaller than our 
results; this is not surprising as the meson resonances included in ref. IQ were the ai(1260), 7r'(1300), a2(1320) and 
a;(1420), the latter three contributing rather little at the free p mass M = irip. On the other hand, the recent kinetic 
theory treatment performed in Ref. p9[ does agree with our findings. However, we would like to stress again that 
our approach consistently accounts for the empirical radiative decays at the same time, which is crucial for reliable 
predictions of low-mass dilepton production to be addressed in the next section. The shift of the pole mass, defined 
by the zero crossing in the real part of the propagator, turns out to be negligible, moving from M = 773 MeV in 
vacuum to M = 776 MeV at T = 150 MeV. 
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FIG. 3. Imaginary part of the /o-propagator (spectral function) in the vacuum (dotted curve), and in a thermal gas including 
the full in- medium self-energies, Eq. (M), for fixed three-momentum q = 0.3 GeV at temperatures T = 120 MeV (long-dashed 
curve), T = 150 MeV (dashed curve) and T = 180 MeV (dotted curve). 



V. DILEPTON PRODUCTION 



The differential dilepton production rate per unit four- volume and four-momentum in hot matter can be decomposed 



as 



dN) 



i+i- 



d^xSq 



L^MH^'iq) 



(35) 



For definiteness, we will focus on e~''e^ pairs in the following. Then the electron/positron rest masses can be neglected 
as compared to their three- momenta, Wgi ^ |p±|, and, to lowest order in the electromagnetic coupling a, the lepton 
tensor takes the form 



L^iy{q) = - 



37r2A/2 



f _ 'Mil) 



(36) 



with the total pair four-momentum q = p^ +p_. Here, we are interested in dilepton production from p decays (or, 
equivalaently, tttt annihilation). Within the VDM, which we have already employed in the evaluation of radiative 
decay widths in sect. O, the corresponding hadronic tensor is directly related to the imaginary part of the retarded p 
propagator in hot and dense matter: 



H^^qo,q';pB,T) = -riqo;T) 






ImD^/{qo,q;fiB,T). 



Using the decomposition Eq. (p4|), the dilepton rate can then be written as 
dN _ ^2(^0)4 fP^q^-T) 1 



d'^xd'^q 



^ 



M2 



- [imD^^ {qo,q;T) + 2 ImDj {qo,q;T)] 



(37) 



(38) 



with the longitudinal and transverse spectral functions given by Eq. (63). 

In Fig. H we display the individual medium effects in the three-momentum integrated rates 

^^ (M^^ fd'qdN 

:{M;T) = / ,.u iqo,Q) ■ 



(39) 



d^xdM'^ V ' ^ y 2qo dHd^q 
for electron-positron production at a fixed temperature of r=150 MeV. The two mechanisms that dominate the net 
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FIG. 4. Three-momentum integrated dilepton production rates at temperature of T = 150 MeV including various individual 



medium effects in the p propagator (see Eqs. (E 



|), (p4|)); dotted curve: free tttt annihilation (i.e., using Ep' 



short-dashed curve: tttt annihilation including finite temperature effects through Bose enhancement factors, Eq. (BU) (using 
Ep '"^ = Ep^^); long-dashed curve: free tttt including pp — > /i resonance formation (using Ep'"^ = Ep.^..^ 4- E ' , ); dashed-dotted 
curve: free tttt including p-K -^ ai resonance formation (using Ep'^ = Ep.^^ -I- Ep^^^); solid curve: free tttt including p-K ^> uj 
resonance formation (using Ep'"^ = Ep.^^ -I- Ep^^^). 
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low-mass dilepton rate are associated with pir ^ u) formation as well as the Bose-Einstein (BE) enhancement of the 
p ^ TTTT decay width. E.g., at M=0.4 GeV, both the w (solid line) and the finite-T (BE) corrections in the pion 
cloud enhance the free tttt annihilation rate by (80-90)% each. Much smaller effects are due to the ai(1260) (30% 
enhancement, dashed-dotted line) as well as the Ki and hi (20% and 15% enhancement, respectively, not shown in 
Fig. 4), the /i being practically negligible. 

In Fig. H we compare the total rate calculated from an incoherent sum of meson reactions and decays [^ (dashed- 
dotted line) with our full result (solid line). Focusing again at M = 0.4 GeV, the latter is increased over the free tttt 
result (=GL, dotted line) by a factor of 3.5, compared to ~ 2 in the Gale-Lichard calculation. The major difference 
arises from the inclusion of the BE-enhancement (and, to a lesser extent, the ai/hi resonances) in the present 
treatment. Also shown in Fig. g| are the results based on the finite temperature part of the calculation by Rapp, 
Chanfray and Wambach (=:RCW, short-dashed curve) ll^, which included oi, Ki resonances and the Bose-Einstein 
enhancement in SpTr^r (note, however, that the hadronic form factors used in Ref. llq] were substantially harder, with 
cutoffs A^paj — h-KpKi = 2 GeV; this entails an overestimation of the radiative ai -^ jn decay width by a factor of 2, 



as elaborated in sect. III). Consequently, in the vicinity of the free p mass, where the major broadening effect is due 
to ai formation and BE-enhancement, the RCW results differ very little from the present ones. Below, w formation 
in piT scattering is responsible for a substantial increase of the emission rate (full curve compared to short-dashed 
curve). 

Another point is finally noteworthy of attention: if one looks at the total spectrum in the vicinity of the r/io-peak, 
one sees that the net signal in our many-body spectral function approach is reduced by about 40% as compared to 
free tttt annihilation. As discussed in the previous section, the shift of the real part is very small (a few MeV; the 
apparent larger shift in the dilepton spectrum is caused by the overall Bose factor in the rate expression, Eq. (|3q), 
which strongly increases towards smaller invariant masses), but the broadening of ~ 80 MeV is appreciable on the 
scale of the free p width of ~ 150 MeV. The resulting peak smearing is a distinct feature of our many-body formalism, 
i.e., the resummation encoded in the denominator in Eq. (B3), which is neither present in the rate calculations of 
Ref. 12^ nor in low-density expansions as performed, e.g., in Ref. [g4|. This feature will be tested against experiment 
soon, as a result of improving the mass resolution in the CERES experiment at CERN-SpS energies, as well as at 
RHIC energies in the PHENIX detector, where the design value for AM/M is at the 1% level. 
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FIG. 5. Comparison of our total, three-momentum integrated thermal dilepton production rate (solid curve) with that 
obtained from a sum of meson decays and reactions according to Ref. El| (dashed-dotted curve) . The dotted curve represents 
the rate obtained from free tttt annihilation, whereas the dashed curve corresponds to the finite temperature part of the 
many-boby calculation of Ref. ||l5|. All curves are for a fixed temperature T = 150 MeV. 
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VI. SUMMARY AND CONCLUSIONS 

Based on a finite temperature (Matsubara) formalism we have highlighted different interaction channels that could 
lead to modifications of the rho properties in a thermal gas of mesons with zero chemical potentials. Using phe- 
nomenological interaction vertices compatible with gauge invariance and chiral symmetry, and including hadronic 
form factors to simulate finite size effects, the free parameters could be tuned to reproduce the empirical hadronic 
and radiative decay branchings rather well. The resulting in-medium p self-energy induces a moderate broadening of 



the p spectral function, somewhat higher than what has been found before in Refs. [£2|,y_5 18 1, but consistent with 
a recent kinetic theory analysis [|l9[. On the other hand, cancellations in the real part inhibited significant changes 
of the in-medium p pole-mass. Corresponding dilepton spectra exhibit a 40% depletion of the p-peak together with 
an appreciable enhancement of a factor of ^ 3.5 below, largely driven by sub-threshold Tip —^ to formation and a 
Bose enhancement in the p ^ tttt width. Our results have now to be combined with those obtained with the baryons 
present, and a time-evolution approach will be coupled with our rates to produce yields that can be compared with 
experiment. 

Finally, for the sake of comparison, we have overlayed our results with the those obtained from an (incoherent) sum 
of meson reactions and decays pl|. Some differences emerge: the earlier calculations in Ref. |2ll do not include the 
contribution that can be associated with the radiative decay of the ai(1260) and /ii(1170); also, the VDM form factor 
used there was temperature- independent. The methods employed in this work are tantamount to the generation 
of a temperature-dependent form factor. Besides that, it seems that the self-consistent many-body treatment does 
not induce large deviations as compared with the results of (incoherent) rate calculations as far as the low-mass 
enhancement is concerned. There are, of course, significant differences around the free p mass, in that the many-body 
calculations lead to a reduction of the p-peak. In fact, these features can be understood in a rather transparent way 
as follows. Schematically, the p spectral function (which directly governs the dilepton rate) can be written in terms 
of the self-energy as 

^"'^^ " |Af2-m2|2 + |ImS,|2 ' ^^°^ 

where we have absorbed the real part of the self-energy in the (physical) p mass rUp. In the low-mass region, where 
rup ^ M and rUp ^ jlmSpj, the denominator is dominated by nip so that 

ImDpiM < nip) oc ^^^ . (41) 

Since ImEp corresponds to a summation of scattering amplitudes times (pion-) density, one immediately recognizes 
the close resemblance to kinetic theory or low-density expansion approaches. On the other hand, in the vicinity of 
the p-peak, where M ~ nip, the denominator in Eq. (HG) is dominated by ImSp so that 

Imi?p(Af ~mp)(x^-^ , (42) 

demonstrating that the consequence of an increase in density is a suppression of the maximum, which cannot be 
straightforwardly casted in a low-density expansion. 

Overall, however, the finite temperature effects found here are still to be regarded as rather moderate. On the 
contrary, the nuclear environment has a much stronger impact on the in-medium rho-modifications at comparable 
densities [ p^|jl^ , p^[ -p^ : e.g., at a temperature T — 160 MeV, where the thermal pion density, n^ = 0.16 fm~^, 
equals normal nuclear density, the p spectral function broadens by about r™'^''(T — 160MeV) ~ 100 MeV, which is 
considerably less than in nuclear matter (r™'^''(p7v = Po) > 300 MeV, as extracted, e.g., from Refs. |lq,pfl]). From 
a theoretical point of view this state of affairs may appear puzzling, since towards the chiral restoration transition, 
which is to be expected around the temperatures considered here, a substantial reshaping of the vector and axial- 
vector spectral distribution must occur: for T > T^, they have to become degenerate as an unavoidable consequence 
of the (approximate) chiral symmetry in the strong interactions. The reason for this discrepancy may be related to the 
fact that in the heat bath the p modifications as calculated in the present article are still hindered by the Goldstone 
nature of the pions. This observation does not apply to p-interactions with nucleons in the finite density case. A 
complete understanding still requires further elucidation. 

Resuming one of the motivations of our analysis, we believe that the understanding of complex strongly interacting 
systems that live in regions of density and temperature far removed from equilibrium resides not only in confronting 
theoretical calculations with experimental data, but also in comparing the theories among themselves. 
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TABLES 



R 


I°J'' 


Ttot [MeV] 


ph Decay 


r°, [MeV] 


r^, [MeV] 


cj(782) 


0^1- 


8.43 


PTV 


~5 


0.72 


fti(1170) 


0-1+ 


-360 


pn 


seen 


? 


ai(1260) 


1-1+ 


-400 


p-K 


dominant 


0.64 


_K'i(1270) 


fl+ 


~90 


pK 


- 60 


? 


/i(1285) 


0+1+ 


25 


PP 


<8 


1.65 


7r'(1300) 


1-0- 


~400 


p-K 


seen 


? 



TABLE I. Mesontc Resonances R with masses rrin < 1300 MeV and substantial branching ratios into final states involving 
direct p's (hadronic) or p -like photons (radiative). 



R 


IF(phR) 


GphE^ [GcV-^] 


KhB. [MeV] 


Kh [MeV] 


r°ft [MeV] 


a;(782) 


1 


25.8 


1000 


3.5 


0.72 


/ii(1170) 


1 


11.37 


1000 


300 


0.60 


ai(1260) 


2 


13.27 


1000 


400 


0.66 


^-1(1270) 


2 


9.42 


1000 


60 


0.32 


/i(1285) 


1 


35.7 


800 


3 


1.67 


7r'(1300) 


2 


9.67 


1000 


300 






TABLE IL Results of our fit to the decay properties of p-h induced mesonic resonances R with masses m_R < 1300 MeV (the 
/i(1285) and 7r'(1300) coupling constants are in units of GeV^'^). 
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